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The Higgs boson is a scalar elementary particle pre- 
dicted by the Higgs mechanism which attempts to 
explain the origin of mass of the elementary particles 
within the Standard Model (SM) [5] . Present experimen- 
tal evidence suggests a Higgs-hke state with the mass of 
« 126 GeV [3]. However, low-mass Higgs states with 
masses of 0(10 GeV) appear in several extensions to 
the SM [3], such as the Next-to-Minimal Supersymmet- 
ric Standard Model (NMSSM) [5,. The NMSSM adds a 
singlet superfield to the Minimal Supersymmetric Stan- 
dard Model, solving the so-called naturalness problem. 
The NMSSM contains two charged Higgs bosons, three 
neutral CP-even bosons, and two CP-odd bosons. The 
lightest CP-odd state, , could have a mass below the 
hh production threshold, avoiding detection at LEP [3]. 
Such a particle could be produced in radiative T — ^ "fA^ 
decays [6] with a branching fraction as large as 10^'' for 



the narrow states T{nS) (where n < 3), depending on 
the A^ mass and couplings making it accessible at 
B-Factorics. Thus, constraining the low-mass NMSSM 
Higgs sector is important for understanding the recent 
LHC discovery j7|. 

Searches for A" decays into il^ n~ [8], r+r" [9], invis- 
ible [TU], and hadronic |TT] final states have been per- 
formed by BaBAR, so far with null results. In particular, 
limits on the product of branching fractions B{T{?,S) — > 
7^°) X B{A° T+T-) have been set at (1.5-16) x 10"^ 
in the mass range 4.03 < m^o < 10.10 GeV [S]. The 
CLEO Collaboration has set limits on the branching ra- 
tio product B{T{IS) -> 7A") x B{A° ^ t+t-) in the 
range 10~^ — 10"* for masses to^o < 9.2 GeV 12 . 

This paper describes a search for decays of the T{1S) 
resonance into a photon and a light scalar CP-odd Higgs 
boson A^ , which then decays into a pair of tau leptons. 
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The T{1S) resonance is produced from the T{2S) reso- 
nance with the emission of two charged pions. The re- 
action chain is e+e" ^ r{2S) Tr+TT-r{lS), r{lS) 
jA°, A° -J> T+T-. We identify the r{lS) by the dipion 
transition; the production and decay of the Higgs candi- 
dates are identified by the photon and the two charged 
tracks from one-prong decays of the two tan leptons. 

This analysis is based on a sample of (98.3 ±0.9) x 10^ 
T{2S) decays collected at the T{2S) resonance with the 
BaBAR detector [T3] at the PEP-II asymmetric-energy 
e+e~ collider at the SLAC National Accelerator Labo- 
ratory. This sample corresponds to an integrated lumi- 
nosity of 14 fb~^. We also use a sample of 28 fb~^ taken 
at the T(3S) resonance for studies of the QED (contin- 
uum) backgrounds and the optimization of the selection 
of the T{2S) dipion transition candidates. T{3S) de- 
cays are rejected by the analysis selection criteria due 
to their different kinematics distributions compared with 
the T{2S) decays, and therefore the T{3S) events form 
a pure high-statistics continuum QED sample. An ad- 
ditional dataset of 1.4 fb~^ taken at a center-of-mass 
(CM) energy 30 MeV below the T{2S) mass is used 
for studies of systematic effects. We use Monte Carlo 
(MC) simulated samples of signal and T{1S) background 
events [HI [15] to tune the selection of the Higgs events. 
The tau-lepton branching fractions are fixed to the val- 
ues of Ref. [in]. The BaBAR detector, including the In- 
strumented Flux Return, the electromagnetic calorime- 
ter, and the tracking and particle identification (PID) 
systems, is described in detail elsewhere [131 [T7]. 

A signal candidate consists of a photon plus four 
charged tracks: tt+tt" from the T{2S) — )■ 7r"'"7r~T(lS') 
transition, and the one-prong decays of the two tau lep- 
tons. The event may contain as many as 19 additional 
photons with laboratory energy greater than 30 MeV, 
mostly from beam-induced backgrounds, but no addi- 
tional charged tracks. Additional signal candidates may 
be formed using these extra photons, but a single final 
candidate per event is selected, as described below. 

We select events where at least one tau lepton decays 
leptonically, resulting in five different combinations of 
tau lepton daughters: ee, e/i, ctt, /i/i, //tt. Events in which 
both tau leptons decay hadronically suffer from signifi- 
cantly larger and poorly modeled backgrounds than the 
leptonic channels, and are therefore excluded. The tau 
lepton daughters are identified using multivariate dis- 
criminants based on the information from all subdetec- 
tors. Typical PID efficiencies are 98% (e), 90% (^), and 
97% (tt) , while the typical pion misidentification rates are 
less than 0.5% (e) and 5% (/u). Requirements on the elec- 
tromagnetic shower shapes of the primary photon candi- 
dates are also imposed to improve the signal purity, and 
events with tt'^ candidates, formed from pairs of photons 
with invariant mass satisfying 100 < m^j < 160 MeV 
and laboratory energy above 200 MeV, are discarded. 

In order to achieve a balanced selection efficiency that 



does not depend strongly on the reconstructed Higgs 
mass niAo (or photon energy), we optimize the selection 
in two Higgs mass intervals: 3.6 < ttiao < 8.0 GeV (L 
range) and 8.0 < m^o < 9.2 GeV (H range). The choice 
of the mass ranges is motivated by the rapidly varying 
signal-to-background ratio at low photon energies (which 
correspond to higher m^o for the signal) , and differences 
in kinematics for each photon energy range. 

The masses of the T{1S) and A" candidates are calcu- 
lated from two primary kinematic variables: 

"^?ccoii = Af2(25^ + ml^ - 2Mr(2S)E^^' , (1) 
m\ = (Pe+e- - - P^f ■ (2) 

Here mrocoii is the recoil mass of the dipion system (which 
peaks at the value of the T{1S) mass for signal), is 
the mass recoiling against the signal photon in the T{1S) 
frame (which peaks near the square of the expected Higgs 
mass, TO^o, and is linear in the photon energy), and P 
denotes the four-momentum. 

In order to reject backgrounds, we train two multilayer 
perceptron neural networks (NN) [TS]: a pion discrimi- 
nant (A/'tt), which describes the kinematics of the process 
T{2S) — > 7r+7r^T(15'), and a tau discriminant (A/^) de- 
scribing the transition T{1S) — ?> 7^°, — > r+r^. Each 
NN uses kinematic variables only weakly correlated to 
m\ or nirocoii- The two discriminants are uncorrelated. 
The pion discriminant M-^ combines nine kinematic vari- 
ables specific to the dipion system [TOl [19] . The discrim- 
inant Mr is a combination of 14 variables: the missing 
energy and the cosine of the polar angle of the missing 
momentum in the event; the extra calorimeter energy in 
the lab frame and the energy of the second most ener- 
getic photon in the CM frame; the net transverse momen- 
tum of the reconstructed signal candidate particles; the 
acoplanarity of the photon relative to the plane formed 
by the two tau decay prongs; the momentum and polar 
angle in the CM frame of the most energetic tau decay 
prong; the invariant mass, vertex probability, and the 
distance between the vertex and the e^e~ interaction re- 
gion of the two tau decay prongs; the angle between the 
signal photon and the most energetic tau decay prong in 
the event. The discriminants are trained using signal MC 
events Y{IS) 7^" in the range 4.0 < m^o < 9.2 GeV. 
The background samples for training are taken from the 
continuum sample for the pion discriminant and from the 
simulated generic T{\S) decays for the tau discriminant. 

Each NN outputs a value N close to -fl for signal and 
to —1 for background. Based on the NN outputs, the 
selection criteria for A/'tt and A/V are chosen to optimize 
e/ {l.b+^/B) [5D], where e is the signal efficiency, and B is 
the expected background yield. We accept signal candi- 
dates if A/V is above a threshold value chosen individually 
for each final state and mass range. The thresholds for 
A/'tt are the same for all final states, but different for the 
two mass ranges. The typical signal efficiency and back- 
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ground rejection estimated from the corresponding MC 
samples are listed in Table |lj 



TABLE I: Typical selection efficiencies (SE) and background 
rejection factors (BR) for the two neural network discrimi- 
nants N-n and Mt in the L (3.6 < m^o < 8.0 GeV) and H 
(8.0 < m^o < 9.2 GeV) ranges. The SE and BR factors are 
relative to the preselection that requires 4 tracks and a real 
photon in the final state, and are averaged over each mass 
range. 



Mass Range 


Discriminant 


SE (%) 


BR (%) 


L 


A/; 


76 


99 


H 


AT, 


72 


99 


L 




80 


97 


H 


ATr 


30 


99 




Due to reconstruction ambiguities, in particular ex- 
tra photon candidates in the event and a large /i-as- 
TT misidentification rate, a fraction of signal and back- 
ground events have more than a single reconstructed can- 
didate. The multiplicity of candidates per event is on 
average 1.8 for the simulated signal samples, 1.6 for the 
generic T(2S) decays, 1.3 for the continuum sample, and 
1.5 for the data. We select a single candidate based on (1) 
the highest value of A/'^r , then, if multiple candidates still 
remain, (2) the highest value oi Afr, and finally, (3) the 
tau decay final state with the highest signal/background 
ratio. 

We further suppress the continuum background by ap- 
plying a cut on the mass recoiling against the dipion sys- 
tem TOrecoii: 

I'Tlrocoil - (TOrGcoil)l < 10 McV , (3) 

where (miocoii) is the expected location of the T{IS) 
peak, determined by the mass difference between the 
T{2S) and T{1S) mesons [21]. The final signal selec- 
tion efficiency varies between 1% and 4.5% (Fig.[l]), and 
is lowest at the highest masses (lowest photon energy). 

We extract the yield of signal events as a function of 
niAo in the interval 3.6 < m^o < 9.2 GeV by performing 
a series of maximum likelihood fits in steps of m^o . We 
perform one-dimensional unbinned extended maximum 
likelihood (ML) fits to the distribution of in the in- 
tervals 12 < mj^ < 72 GeV^ (L range) and 49 < < 89 
GeV^ (H range). The fit intervals overlap to provide suf- 
ficient sidebands for all values of m^o. The likelihood 
contains contributions from signal, which is expected to 
peak near the Higgs mass squared, and from the smooth 
background function, arising from continuum and radia- 
tive leptonic T{1S) backgrounds. We search for the A'^ in 
varying mass steps that correspond to approximately half 
of the expected resolution on m^o, as described below. 
A total of 201 mass points are sampled. 

We use signal MC samples T{2S) -> 
7r+7r~T(lS'), T{1S) — generated at nine (seven) 



FIG. 1: Signal efficiency as a function of m^o in the L and H 
mass ranges (separated by the vertical blue line). 



different values of uiao in the L (H) mass range to de- 
termine the signal probability density functions (PDFs) 
in and selection efficiencies. We interpolate these 
distributions and efficiencies between fixed m^o points, 
correcting for known small differences between data and 
MC simulations. The signal PDFs only include events 
in which the simulated signal photon, tt+tt" tracks, 
and the tracks from the decay of the tau leptons are 
correctly reconstructed. The signal efficiency, however, 
includes contributions from events in which one of the 
final-state particles may be misreconstructed. The signal 
PDF is modeled as a Crystal Ball function J22]. We 
parameterize the background PDFs as 

/(z) = (Erf [a{z - z,„i„)] + 1) ^3(2) (L range) ; 
f{z) = (1-—) exp[C2(z)] (H range) , 

where z = m?x, Cn{z) is an nth-order Chebyshev poly- 
nomial (with different parameters for each mass range), 
a and /3 are threshold parameters, and Zmin = 12.5 GeV^ 
and Zmax — 88.9 GeV^ are determined by the kine- 
matic end points of the photon energy spectrum in the 
luminosity-weighted mixture of simulated generic T(\S) 
decays and continuum sample. A common background 
PDF describes all five decay channels adequately well. 
Parameters of the Chebyshev polynomials Cn{z) and the 
threshold parameters a and (5 are determined from a fit 
to the data distribution of with the signal yield fixed 
to zero. This accounts for uncertainties in the modeling 
of the radiative T{1S) decays and additional backgrounds 
that may not be well described by the continuum sample, 
such as two-photon events with low photon energy (high 
m\). 

For each m^o hypothesis, we determine two param- 
eters: the background yield A^bkg and the signal yield 
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A^sig- In the H range, we also allow the two coefficients 
of C2{z) and the parameter P to vary in the fit. The 
fit is performed simultaneously over the distributions in 
each T+r~ decay channel, taking advantage of the dif- 
ference in the signal-to-background ratios over the decay 
channels. The fraction of events in each channel is fixed 
from MC samples for signal, and from the luminosity- 
weighted mixture of simulated generic T(1S) decays and 
continuum sample for the background. 

Each fixed nuisance parameter in the fit is varied ac- 
cording to its uncertainty; correlations between parame- 
ters are taken into account. The systematic uncertainties 
for this measurement can be divided into two categories: 

• Additive errors: uncertainties on the event yield, 
which do not scale with the number of recon- 
structed signal events. These include uncertainties 
of the parameters fixed in the fit (PDF shape pa- 
rameters for signal and backgrounds) and a small 
bias in the ML fit. These uncertainties reduce the 
significance of any observed signal |23j . 

• Multiplicative errors: uncertainties that scale with 
the number of reconstructed signal events. These 
include uncertainties on the reconstruction effi- 
ciency, the ML fit bias which scales with the true 
number of signal events, the uncertainty in the 
number of produced T{2S) mesons, and the un- 
certainty in the branching fraction of T{2S) — > 
7r+7r-r(lS'). 

We compute the average bias of the ML fit for a set 
of generated N^ig values using a large ensemble of sim- 
ulated pseudo-experiments. In each pseudo-experiment, 
the signal events are fully simulated, and the background 
events are sampled from their PDFs. We determine the 
fit bias that is independent of TVgig and is part of the 
additive uncertainties, as well as the bias which scales 
linearly with iVgigj can be thought of as a "fit in- 
efficiency" , i.e. a relative correction to the signal recon- 
struction efficiency. The bias arises from imperfections in 
modeling of the signal PDFs, from events in which signal 
candidates are misreconstructed, and from low-statistics 
properties of the ML estimators. We see that a correction 
of 3.1 ± 1.1% (L range) and 7.6 ± 1.4% (H range) has to 
be applied. The additive parts of the fit bias are small. 

The signal efficiencies determined in MC simulations 
are corrected by several multiplicative effects: 

• Tracking and dipion selection efficiency. These 
corrections and their uncertainties have been 
determined [10] using a clean sample of of 
four-track final states from decays T{2S) — > 
7r+7r-r(15'), T{IS) The data/MC ratio 
of 0.97 ± 0.02 includes the uncertainties due to the 
number of produced T{2S) events, dipion branch- 
ing ratio T{2S) — )• 7r"'"7r~T(lS'), dipion reconstruc- 
tion efficiency, the efficiency of reconstructing two 



additional (energetic) charged tracks, trigger un- 
certainties, and the selection efficiency for the pion 
discriminant M-n . The uncertainty is dominated by 
the error on T{1S) l^~^IJ-~ branching ratio [TB]. 



• Photon selection efficiency. A correction of 0.967± 
0.017 is determined from a high-statistics e'^e~ — >■ 
77 sample in which one of the photons converts in 
the inner detector material to produce a detectable 
e+e^ pair ilOj. 

• Neural Network selection efficiency. We evalu- 
ate the systematic uncertainty due to possible 
data/MC differences in the distributions of the 
NN discriminant Afr using an inclusive background 
sample. We select signal-like events that pass the 
requirement Afr > and compute the ratio of par- 
tial selection efficiencies for the actual Afr thresh- 
olds for the data and the background. These ratios 
are 1.038 ± 0.013 (L range) and 0.991 ± 0.014 (H 
range) . 

The total correction to the efficiency is a product of all 
corrections discussed above: 

edata/£MC = 0.943 ± 0.031 (L range) ; 
edata/£MC = 0.859 ± 0.033 (H range) . 

We compute the statistical significance of a particular 
fit centered at m^o asS = ^2 log(£niax/'Co), where £max 
is the maximum likelihood value for a fit with a free signal 
yield, and £o is the value of the likelihood for A^sig = 
0. Including additive systematic uncertainties, the most 
significant upward fluctuations occur at m^o = 6.36 GeV 
with S = 2.7(7 (Fig. [2}i) and m^o = 8.93 GeV with S ^ 
3.0(7 (Fig. [2]3). Fluctuations of -1-3(7 or higher occur in 
7.5% of pseudo-experiments that simulate a scan of 201 
mass points with an average correlation of 94.5% between 
adjacent points, as observed in our dataset. Therefore, 
we conclude that no significant signal is found. 

Since we do not observe a significant excess of events 
above the background, we set 90% confidence level (C.L.) 
Bayesian upper limits on the product S(T(15') 7^°) x 
Bi^JP —7- T+T^), computed with a uniform prior for 
A'sig > 0. The limits are shown in Fig. |3 Systematic 
uncertainties on iVgig and Sdata are included assuming 
their likelihood profiles are Gaussian [53]. 

We combine our results with the previous limits on the 
branching ratios i3(T(3S') ^ 7^") x B{A^ r+r") [5] 
to set a limit on the Yukawa coupling of the 6-quark 
to the A". The branching fractions B{T{nS) jA'^) are 
related to gb through [6] [241125]: 



B{r{nS) 7^°) 
BirinS) ^ 1+1-) " y^Tta 



glGprnl 



QCD 



1 - 



(4) 
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FIG. 2: Fits to the distributions in (a) L and (b) H 
ranges for the two particular m^o points that return the 
largest upward fluctuations: (a) m^o = 6.36 GeV and (b) 
m^o — 8.93 GeV. The red solid line shows the signal PDF, 
the green dot-dashed line is the background contribution, and 
the blue solid line shows the total PDF. The top plots show 
the fit residuals normalized by the error (pulls). The signal 
peaks corresponds to a statistical significance of (a) 2.7a and 
(b) 3.0cr. 
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FIG. 3: 90% C.L. upper limits for B{T{1S) ■yA'-'')xB{A° 

T + T-). 



as a function of my^o are shown in Fig. [4j They rule 
out much of the parameter space preferred by NMSSM 
gt — tan /3 cos 6*^1 > 1, where tan/3 is the ratio of the 
Higgs vacuum expectation values and cos 9a is the frac- 
tion of the non-singlet component in the CP-odd A" [5]. 
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where / = e or /i, a is the fine structure constant com- 
puted at the scale m'f(nS): Gf '^s the Fermi constant, and 
J'QCD includes the to^o -dependent QCD and relativistic 
corrections to B{T{nS) — 7^°) [25] and the leptonic 
width of T{nS) [55]. To first order in as, the corrections 
range from to 30% [25^ and may have large uncertain- 
ties [37]. 

We combine our data on T{1S) ^A^ with the BaR4R 
results of Ref. [9^ using the full likelihood functions for 
at each m^o point from this analysis, and a Gaussian ap- 
proximation for the g^ likelihood from Ref. . The com- 
bined upper limits on the product g^ x B{A'^ — t+t~) 



FIG. 4: 90% C.L. upper limits for Yukawa coupling x 
B{A''^ — >■ T^r~). Shown are combined BABAR results (red 
solid line), results from this analysis only (dashed green line), 
the previous BABAR measurement [5] (dotted blue line), and 
results from the CLEO experiment [12] (dot-dashed black 
line). The shaded vertical bar shows the region around X6(2P) 
mass excluded from Ref. |5]. 



In summary, we find no evidence for the radiative de- 
cays T{1S) — >■ in which A^ decays into a pair of tau 
leptons, and we set 90% C.L. upper limits on the prod- 
uct of branching fractions B(T{1S) — 7^°) x B{A" — 
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T+T") in the range (0.9 - 13) x 10"^ for 3.6 < tuao < 
9.2 GeV. We also set 90% C.L. upper limits on the prod- 
uct gl X B{A" T+T-) in the range 0.09 - 1.9 for 
m^o < 9.2 GeV. Our limits place significant constraints 
on NMSSM parameter space. 
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